Charged Kaon Mass Measurement using the Cherenkov Effect 
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Abstract 



^ The two most recent and precise measurements of the charged kaon mass use X-rays from kaonic atoms and report uncertainties of 
14 ppm and 22 ppm yet differ from each other by 122 ppm. We describe the possibihty of an independent mass measurement using 
the measurement of Cherenkov light from a narrow-band beam of kaons, pions, and protons. This technique was demonstrated 
using data taken opportunistically by the Main Injector Particle Production experiment at Fermi National Accelerator Laboratory 
which recorded beams of protons, kaons, and pions ranging in momentum from -|-37 GeV/c to -|-63 GeV/c. The measured value 
is 491.3 ± 1.7 McV/c^, which is within \A<j of the world average. An improvement of two orders of magnitude in precision would 
make this technique useful for resolving the ambiguity in the X-ray data and may be achievable in a dedicated experiment. 



1. Introduction 

The charged kaon mass is an important input in deter- 
mining the CKM matrix element Vus from measurements 
of the branching ratio of K'^ — >■ Tr^e+i^. The value of the 
charged kaon mass reported by the Particle Data Group is 
493.677 MeV/c^ with an uncertainty of 26 parts per million 
(ppm) [1] . This value is a weighted average of six measure- 
ments but is dominated by the two most recent and precise 
measurements from Denisov [2] and Gall [3] which measure 



X-ray energies from kaonic atoms. While these measure- 
ments report uncertainties of 14 and 22 ppm they differ 
by 122 ppm (4.6cr). In this article we explore one possibil- 
ity to resolve this discrepancy using an independent tech- 
nique for measuring the charged kaon mass based on the 
Cherenkov effect. The well known pion and proton masses 
are used as references. The technique is demonstrated us- 
ing data taken opportunistically using the Ring Imaging 
Cherenkov (RICH) sub-detector of the Main Injector Par- 
ticle Production (MIPP) experiment at Fermilab [4]. 



2. Measurement Concept 

Cherenkov light is emitted when a relativistic charged 
particle of mass m, momentum p, and speed /3 = 
1 / y^l + (m/p)"^ travels through a radiator volume of in- 
dex of refraction n with /3 > 1/n. (As is customary in 
high energy physics we work in units in which the speed 
of light, c, is 1.) Neglecting dispersive effects for the mo- 
ment, the light is emitted in a cone at an angle 6 given by 
cos^ = 1/Pn [5,6] which is approximately 



the differences in the Cherenkov angles are largest, while 
staying above proton threshold. 

In a RICH detector, the angle 6 can be determined on 
a track-by-track basis from the pattern of Cherenkov pho- 
tons recorded. However, the light for a single ring will be 
distributed about the central angle 6 due to the variation 
of the index of refraction of the radiator medium over the 
wavelengths at which Cherenkov photons are produced. 
This gives as contribution to the uncertainty in the average 
angle 6 determined from a single track of 



(1) 



for small angles. Now, consider two particles i and j with 
identical momenta p but different masses mj, rrij and speeds 
(3i, Pj. They will emit Cherenkov light at angles 6i, 9j which 
are related by the expression 



Mi - 0jO] = 2(A - Pj). 

In the relativistic limit p ^ m, 
combined with Eqn. 2, gives 
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1 - This, when 
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where we have neglected the small difference between 39 
and 9. If the particles i and j are pious, protons, and kaons, 
we have two independent angle-squared differences that can 
be measured 
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Using these, the kaon mass is given by 
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where Ap^ = m^ — mj. Notice that for a monochromatic 
beam in the absence of dispersion the index of refraction 
n and momenta p drop out. The kaon mass can be deter- 
mined, in principle, through measurements of the pion and 
proton masses and the Cherenkov angles of the three par- 
ticles. The proton and pion masses are known to 0.9 ppm 
and 2.5 ppm respectively and will not be the limiting fac- 
tors in the experiment. 

Using Eqn. 5 we can estimate the uncertainty in m\ 
measured using this method as: 



+ 



2 

K \ J2 , 
CTa + 



4/ 



92 _ Q2 

A2 a 2 

2 pK 2 ,a2 2 I q2 ^K-k 2 



A2 



(6) 



The first two terms are due to the uncertainties in the pion 
and proton masses and are small. The third term grows 
with momentum and suggests that it is best to conduct 
the measurement at as low a momentum as possible where 
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whcire (5„ is the amount of dispersion over the photomul- 
tiplier tube (PMT) wavelength acceptance and Nh is the 
number of PMT hits in the ring. 

In a beamline, particles will be accepted if their momen- 
tum lies in a narrow window about some central value p. 
The finite size of this momentum acceptance window intro- 
duces an additional uncertainty in the average angle mea- 
sured from a single track. Averaging over rings, the mo- 
mentum spread contributes an uncertainty 
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to the measurement of 9, where Op is the spread in the 
beam particle momenta about their central value. We take 
the specific case of an experiment using CO2 as the radi- 
ator (n = 1.00045, (5„ = 3 X 10~^) and a beam of central 
momentum 40 GeV/c with a width of ijp/p=0.01. Under 
these assumptions, ae^ values are in good agreement with 
observed widths of ring radius distributions. This gives a 
statistical uncertainty in the kaon mass of 
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where the first term results from uncertainties in the pion 
mass, the second from uncertainties in the proton mass, and 
the final three terms result from uncertainties in the angles 
9„, 9k, and 9p with A^, Nk, arid Np being the number 
of millions of pion, kaon, and proton rings recorded. Using 
Eqn. 9 we find that the uncertainty is minimized if 32% of 
the data is collected using protons, 23% using pions, and 
45% using kaons. This result is only weakly dependent on 
momentum as shown in Figure 1 which plots the expected 
statistical precision of the charged kaon mass as a function 
of momentum choice and total number of rings recorded. 
With 10 million rings at 40 GeV/ c, we expect a statistical 
precision of 30 ppm using this technique. 

3. RICH detector overview 

The RICH detector used by MIFF was built by the SE- 
LEX Collaboration [7] for use in that experiment. We sum- 
marize here only the most important features of the de- 
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Fig. 1. Expected statistical uncertainty for kaon mass measurement. 

tector as deployed for the MIPP experiment and refer the 
reader to [8-11] for details. 
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Fig. 2. Schematic of the RICH detector [8]. 

The geometry of the RICH counter is shown in Figure 2. 
The detector was constructed from a low carbon cylindrical 
steel vessel 10.22 m in length and 93 in. in diameter with 
a wall thickness of | in. The ends were sealed with 1.5 in. 
thick aluminum flanges that were cut out to hold thin beam 
windows at each end and a photomultiplier tube holder 
plate at the upstream end. The vessel was wrapped by a 
water line carrying chilled water and 15 cm of building 
insulation to regulate the temperature. 

A2.4mxl.2m array of 16 mirrors mounted at the down- 
stream end of the counter focused Cherenkov light on an ar- 
ray of 1/2 in. photomultiplier tubes. On average the mirrors 
had a radius of curvature of 1980 cm with variations less 
than 5 cm and a reflectivity of more than 85% at 160 nm. 

MIPP used CO2 gas held at just above atmospheric pres- 
sure as the radiator. At STP, CO2 has an index of refrac- 
tion of 1.00045 at A = 300 nm giving thresholds oi p — 
4.5 GeV/c for pions, 17 GeV/c for kaons, and 31 GeV/c for 
protons. The gas temperature and pressure were continu- 
ously monitored enabling calibration of the index of refrac- 
tion. A /? = 1 particle produced a ring of radius 29.5 cm and 
an average of 30 PMT hits providing 3ct -k/K separation 
up to 80 GeV/c and 3cr p/K separation up to 120 GeV/c. 
Figure 3 shows sample event displays for each beam species 
showing PMT hits and reconstructed rings. 

The PMTs were mounted in a hexagonal array behind 
a 2 mm thick quartz window which provided a gas-tight 
seal between the phototubes and the radiator volume. They 
were read out in threshold (on/off) mode. Two models 
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Fig. 3. Sample event displays of 40 GeV pion (top), kaon (second 
from top), and proton (second from bottom) rings in the RICH 
counter. Small circles indicate hit PMTs. The large dashed circles 
show the rings reconstructed from the PMT hits. The bottom plot 
is the distribution of reconstructed ring radii for 40 GeV. 

of PMTs were used: Hamamatsu R760 and Russian made 
FEU60. The FEU60 tubes were coated with a wavelength 
shifter to match the acceptance region of the R760's. The 
R760 (FEU60) tubes have a maximum quantum efficiency 
of about 25% (11%) at 350 nm. Of the available 89 PMT 
columns, 68 were used, with R760's installed in 15 columns 
and FEU60's installed in 53 columns. The front-end elec- 
tronics of the RICH detector was re-designed and rebuilt 
by the MIPP collaboration. 



4. Data Analysis 

In outlining the measurement concept, we made several 
simplifying assumptions which must be accounted for dur- 
ing analysis. Due to dispersion, Cherenkov light is not ob- 
served at a single angle, but as a distribution across several 
angles. In our analysis, we measured the average PMT oc- 
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cupancies for pion, proton, and kaon rings as a function of 
Cherenkov angle in data and compared them to calculations 
which incorporated the kaon mass as a free parameter. In 
total 12 million rings were recorded for this measurement 
using positively charged beams of pions, kaons, and protons 
ranging from 37 GeV/c to 63 GeV/c in momentum. 

In the data, wire chambers upstream and downstream 
of the RICH counter were used to reconstruct the particle 
trajectory and predict the center of the RICH ring. Using 
that prediction for the center position, a Cherenkov angle 
was assigned to each PMT and the average probability for 
a PMT to fire was computed as a function of Cherenkov 
angle. This procedure was done separately for each PMT 
type, mirror, and momentum setting. As the measurement 
was made with beam particles, light struck only the two 
central mirrors, labeled 8 and 9. 

To compare to the data, we calculated the expected oc- 
cupancy of each PMT in the RICH array for pion, proton, 
and kaon rings. The calculation starts with the number of 
Cherenkov photons produced per unit path length between 
wavelengths A and X + dX [5,6]: 
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where a is the fine structure constant. These photons 
travel through the radiator medium, reflect off mirrors, 
pass through a quartz window, and are collected by a re- 
flective cone before they are incident on the photodetector. 
The transmission probability for a Cherenkov photon of 
wavelength A for each of these steps is plotted in Figure 4 
along with the photodetector efficiency. Combining these 
factors, the average number of photoelectrons detected by 
PMT i is given by 
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where e(A) is the product of all wavelength-dependent effi- 
ciency factors, /i(A) is the absorption coefficient of CO2, Fl 
is the mirror focal length, and Gi is the geometric accep- 
tance of the i^^ PMT. Scattering of light from angular bin 
6c to bin 9 is accounted for by the function S{9, 9c)- The 
light scattering was modeled as a Gaussian of width a with 
three components. These include an intrinsic term indepen- 
dent of wavelength and a wavelength-dependent dispersive 
term [12]. A third term accounts for multiple scattering of 
the beam particle in the RICH radiator: 
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where ctq is the intrinsic scattering width and (T„ is the 
dispersive scattering width. The calculated shapes are im- 
plicitly functions of the particle masses through /3. Using 
Poisson statistics, the probability for a PMT seeing Np^ 
photoelectrons to be on is: 
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Fig. 4. Overlay of all efficiency functions incorporated into the calcu- 
lation of detector photoelectrons as a function of angle. Also shown 
are the initial Cherenkov photon production spectrum and final pho- 
toelectron yield normalized to peak at 1 for inclusion in the plot. 



P = l-exp(-Arpe-6), 



(13) 



where b accounts for the PMT dark noise rate. 

The measured and calculated PMT occupancies were 
compared using a statistic. In the comparison several pa- 
rameters were allowed to vary to minimize . These were: 

i. Intrinsic detector smearing width, ctq 

ii. Dispersive smearing width, an 

iii. Density ratio scaling factor 

iv. Density ratio offset 

V. Level of air contamination in CO2 

vi. Central pion momentum 

vii. Width of momentum distribution 

viii. Kaon mass 

The density ratio is defined as the gas density inside the 
RICH divided by density at STP and was varied to account 
for uncertainties in the calibration of the pressure and tem- 
perature monitors installed in the counter. 

The momentum acceptance of the beamline was modeled 
as a Gaussian and uncertainties in the mean and width of 
this acceptance were incorporated into the calculation. The 
average proton and kaon momenta accepted by the beam- 
line differ slightly from the average pion momentum due to 
variations in particle production at the secondary target. 
This effect was incorporated using a simulation of the cop- 
per secondary target using the FLUKA [13,14] production 
model. 

After all known effects were accounted for there were 
0(10%) differences between the measured and calculated 
PMT occupancies, due presumably to our incomplete un- 
derstanding of the factors which control Cherenkov photon 
production and transport in our detector. To account for 
these differences, the uncertainties in the PMT occupan- 
cies were increased until ^ jndj = 1 was achieved at the 
minimum. The scaling factor ranged from about 25 to 100 
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for the different combinations of PMT type, mirror, and 
beam momentum. Systematic uncertainties overwhelm the 
statistical uncertainties for this data set. A sample fit for 
40 GeV/c data is shown in Figure 5. Percent difference be- 
tween data and predicted occupancy is shown in Figure 6. 
They agree within 5-15% over the angular rangue illumi- 
nated by kaon rings. This difference is typically 15 times 
larger than the statistical uncertainty in each angular bin. 
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Fig. 5. Comparison of measured PMT occupancies to best-fit pre- 
dictions for 40 GeV data set [15]. Pions are shown in the top row, 
protons in the middle row, and kaons in the bottom row. In each 
row, the left panel shows the results for R760 PMTs while the right 
panel shows the results for FEU60 PMTs. 
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Fig. 6. Percent difference between data and predicted occupancy for 
40 GeV/c kaons. Hamamatsu tubes are on the left, FEU tubes are 
on the right. 



Fig. 7. Summary of kaon mass results for data sets using various 
beam momenta and mirrors. At each momentum setting, results for 
mirror 8 are plotted shifted slightly to the left, while results for 
mirror 9 are plotted shifted slightly to the right. The solid line shows 
the final kaon mass result, and the dotted line is the accepted value 
from the PDG. 

Kaon mass results for each data set are shown in Figure 7. 
Using just the low momentum data sets results in a kaon 
mass measurement of 491.9 ±1.1 MeV; using just the high 
momentum data sets gives 486.7 ± 3.0 MeV. These results 
agree within uncertainties and are combined to give a final 
result for the charged kaon mass of 491.3 ± 1.7 MeV (3500 
ppm) . The convention of scaling error bars when combining 
measurements so that jndf = 1 is followed resulting in 
a larger final uncertainty. This result is within 1.4cr of the 
PDG value [1]. 

See [15] for a detailed discussion of the analysis presented 
in this paper. 

5. Suggestions for Improvement 

In the analysis above we achieved an understanding of 
the PMT occupancies as a function of Cherenkov angle at 
the 10% level resulting in a final uncertainty in the kaon 
mass of 3500 ppm. If the uncertainties in the PMT occu- 
pancies could be reduced to 1% the contribution of this 
source of uncertainty in the kaon mass would be reduced 
to the level of the statistical uncertainties and make this 
technique useful for resolving the 122 ppm discrepancy in 
the X-ray measurements. 

The largest uncertainty in our calculation of the PMT 
occupancy as a function of angle results from knowledge of 
the acceptance of the PMT array which was complicated 



by cross-talk in the readout electronics - a problem which 

manifested itself during the high rate conditions the detec- 
tor was operated under for this measurement. Occasionally, 
a single PMT hit would cause all 16 channels sharing the 
same readout board to register a hit. In our analysis we re- 
jected PMT hits that appeared to be due to cross-talk. As 
the cross-talk was difficult to model, its affect on the PMT 
acceptances has large uncertainties accounting for most of 
the ^10% uncertainty in PMT occupancy as a function of 
angle. 

Additional uncertainties in the occupancies resiilt from 
our treatment of the detector response as a function of 
wavelength. A measurement of the PMT response on a 
tube-by tube basis woiild have required dis-assembly of the 
detector and was not undertaken, rather we used a single 
average response for all PMTs. As the response of individ- 
ual PMTs typically differ from the average by about 15%, 
this introduces an uncertainty of ^6% in the calculation of 
occupancy versus angle. Likewise, we assumed identical re- 
flection efficiencies for the two reflecting mirrors. Taken to- 
gether, the uncertainties in the PMT acceptances and spec- 
tral response of the optical system contribute 3400 ppm to 
our measurement of the kaon mass. 

During the run the temperature and pressure of the ra- 
diator gas varied considerably and the index of refraction 
had to be calibrated several times a day. We estimate that 
uncertainties in the knowledge of index of refraction con- 
tributed 500 ppm in our kaon mass measurement. The 
spread on the beam momentum about its central value con- 
tributed roughly 50 ppm to the kaon mass measurement. 

To achieve 40 ppm precision with this technique each of 
the uncertainties listed above would need to be addressed. 
The largest improvements are to be gained by using a sim- 
pler optical system which could be completely characterized 
and operated at a more consistent pressure and tempera- 
ture free from electronic cross-talk problems. For example, 
the optical system would be simpler if only a single primary 
mirror could be used and the transmission windows elim- 
inated. Use of a fine-grained optical detector would allow 
for the device to be made much more compact allowing for 
better control of pressure and temperature. A more com- 
pact device could be built to run with large excursions in 
gas pressure and temperature allowing for the variation in 
the index of refraction with temperature and pressure to 
be studied in situ. Finally, a charge based readout (rather 
than the on/off readout used in the above analysis) woiild 
allow for better pixel-by-pixel knowledge of the detector 
response. This configuration could be achieved, for exam- 
ple, by replacing the secondary mirror in the MIPP beam 
Cherenkov counter with a GEM-based photo-cathode [16] 
which would afford 30 micron resolution of the Cherenkov 
rings in a device with is 45.7 cm in diameter and 22.9 m 
long. A complete assessment of the performance of such a 
device would require a significant research effort and awaits 
an experimental proposal. We hope that this initial work 
will inspire these future investigations. 

In conclusion, we have used the Cherenkov effect to mea- 



sure the charged kaon mass with 3500 ppm precision us- 
ing an existing RICH detector. We believe that there are 
significant opportunities for a future experiment dedicated 
to the kaon mass measurement to reduce systematic un- 
certainties to the level where this tcchniqiie may be useful 
to resolve the 122 ppm discrepancy in the X-ray measure- 
ments of the charged kaon mass. 
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